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Abstract 

We study the effect of interlayer tunneling in the gauge theory describing a quasi- 
two-dimensional paramagnetic metal close to a second-order or weakly first-order 
antiferromagnetic phase boundary. In that theory, two species of fermions have op- 
posite (rather than equal) charges with respect to the gauge field. We find that 
single-particle interlayer tunneling is suppressed at low energies. The effect of pair 
tunneling is analyzed within the (3 — d) expansion. The resulting phase diagram has 
superconducting and non- Fermi-liquid normal phases, and so is compatible with that 
of the copper-oxide superconductors. 
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Recently, a non-trivial fixed point and the associated non-Fermi-liquid behavior 
have been found for a system of two-dimensional (2d) fermions interacting with a 
transverse gauge field see also [p] — QTD| . The idea of non- Fermi-liquid metal has 



long been proposed [[UJ to apply to the copper-oxide superconductors. In that case, a 



study of interactions between two copies (layers) of such 2d systems is important for at 
least two reasons. First, the real materials are three-dimensional and one would like to 
know if interlayer interactions destroy the fixed point. Second, interlayer tunneling of 
pairs of electrons, in cases when single-electron tunneling is suppressed by non-Fermi- 
liquid effects, was suggested as the mechanism for high-T c superconductivity |T5| ]; 



it seems worthwhile to explore this possibility within the gauge theory framework. 

A priori, one can imagine various types of charge assignments for fermions with re- 
spect to the new gauge field. The correct one, for a real system, should be determined 
from microscopic considerations. In the model considered in refs. @ all fermions 
have the same charge. Here we want to consider a different model, in which there are 
two species of fermions with opposite values of the charge. This charge assignment 
can be derived from microphysics if the antiferromagnetic (AF) correlation length 
extends over at least a few lattice spacings 0, [7]]. Such situation naturally arises 



in systems close to a second-order or weakly first-order AF phase boundary. In the 
copper oxides, it is likely to occur at doping concentrations lower than those for which 
T c is maximal — the "underdoped" materials; however, the model itself may remain 
correct even at the optimal doping concentrations. We do not have anything to add 
to the derivation of this model here, but concern ourselves with its phenomenology, 
namely, the phase diagram. 

The main new element of our theory is that we outfit the model with the oppositely 
charged fermions with local interlayer (and in-layer) interactions. The difference in 
the fate of interlayer interactions between the two types of models discussed above 



is indeed dramatic. In the model of refs.Q-Q, pair tunneling is irrelevant at low 
energies, at least in the leading order of the approximations employed. In our model, 
while single-particle tunneling is suppressed at low energies by the spin gap, pair 
tunneling is relevant and together with the local in-layer interaction determines the 
phase diagram of the system. 

Working in the leading order of the e = 3 — d expansion [|j], we find two phases: 
superconducting and non-Fermi-liquid normal. In the leading order, the anomalous 
dimension of the fermions in the normal phase coincides with that found in refs.Jl], [|; 
we argue that to hold to all orders, despite the presence of a non-vanishing repulsive 
local interaction. A phase transition with temperature between the two phases ap- 
pears possible. Thus, our model, besides being to some extent motivated from the 
microscopic standpoint, has a phase diagram compatible with that of the copper-oxide 
superconductors. 



The new gauge field is supposed to arise as a result of spin-charge separation, 
with electron dissociating into two new quasiparticles. In the Neel phase on the square 
lattice we would have on each site of sublattice A 

^Aa = Z a f A (1) 

and on each site of sublattice B 

IpBa = z a fl , (2) 

where ip is the electron operator and a is the spin index. The smoothly varying 
bosonic fields z and z represent the staggered spin variable 

n = z^crz = —z^az (3) 

and therefore are related as 

~z a = e ap z* (4) 



(up to an inessential phase factor) |I4| , [7j. The redundancy in the number of compo- 
nents of z is the gauge redundancy, apparent in (JlJ) . Because z and z* have opposite 
charges with respect to the gauge field, so do and The fermionic fields Ja and 
Jb also carry the usual electric charge, the same amount for both species. In the AF 
phase, the new gauge symmetry is spontaneously broken by a non-zero expectation 
value of z. 

We are interested, though, in the paramagnetic phase. It is more difficult to 
visualize. However, the charge assignments of the fields cannot change across the 
phase boundary, at least as long as the correlation length stays sufficiently large. On 
both sides of a second-order or weakly first-order phase transition a system should 
be described by the same effective theory, the only difference being in the values of 
the parameters. Therefore, we describe the paramagnetic phase by the same theory 
containing two oppositely charged fermions (the index A,B now merely labels the 
species) and the field z, interacting with the gauge field a M . However, the expectation 
value of z is now zero and all its components are massive ("paramagnetic magnons"). 
The Lagrangian density describing the gauge interactions of these fields in a single 
layer is given in refs.[|l4], [7|; our notations follow ref. 0. By the same logic as above, 
even as we move further away from the phase boundary, the interaction of the low- 
energy fields — the fermions and the gauge fields — should retain its form (see Eq.@ 
below). 

When there are two such layers, fields in each of them have their own gauge 
transformations. The only gauge-invariant interactions between the layers are those 
that correspond to tunneling of objects neutral with respect to the gauge fields. 
Thus, a single fermion can only tunnel accompanied by a quantum of the z field. The 
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corresponding lowest- dimension operators are 

O sA = fl A f2AZ ] 2 Zl (5) 

and similar ones for B fermions or A^B tunneling; 1,2 is the layer index. At frequen- 
cies much lower than the spin gap (which determines the mass of z) and distances 
much larger than the AF correlation length, the tunneling process described by (|5]) 
is suppressed. In that regime, the field z can be integrated out leaving only infrared- 
finite renormalizations of local interactions, such as Eqs.(^j), (0) below. 

Pairs of fermions can tunnel without z quanta. The pair tunneling operators are 

Op = fhfl B f2Bf2A (6) 

and its hermitean conjugate. Within the effective low-energy theory, these are viewed 
as local operators in the 2d space and time, analogous to the reduced BCS Hamilto- 
nian. (For example, in the effective theory there is no objection to taking two holes 
at the same space point, as in (^).) The pair tunneling leads to a renormalization of 
the in-layer four-fermion interactions 

Orl = flAflBflBflA , O r2 = , (7) 

which we then have to include in our analysis. The particle-hole tunneling operators 
not involving z fields are 

fiAfiAfsAfsA, fhfiAfliifw etc. (8) 

but these will turn out to be irrelevant at small values of the corresponding coupling 
constants. 

The physical situation we want to describe is as follows. We assume that the gauge 



field has become deconfining due to effects of the gap less fermions [16]. Because in 



our model A and B fermions have opposite gauge charges, A and B particles moving 
with opposite momenta make parallel currents and are attracted to each other. This 
attraction increases their chances to be near each other and thus effectively enhances 
the local interactions (P), (0). We then have to understand if this enhancement can 
be limited by a repulsive effect of the local interactions themselves. 

There may be different ways to describe the physics outlined above. Here we use 
the perturbation theory built with the resummed gauge propagator || , incorporating 
the Landau damping, but with the usual free fermion propagator (in contrast to the 
1/N method of refs.[|IJ, |3|). The advantage of this approach is that the strongest 
infrared effects of the local four-fermion interactions are easy to isolate. These are 
the infrared-singular renormalizations coming from the BCS channel. Another useful 
device is the (3 — d) expansion ||. At d — > 3, the infrared divergences due to the 
gauge interaction are also logarithmic, and the logarithms can be conveniently handled 
by the renormalization group. Rather than using a modified fermion propagator, 
the (3 — d) expansion introduces an additional renormalizable parameter — the Fermi 
velocity. Such a trade-off is familiar from the A0 4 scalar theory: in the (4 — d) 
expansion A is a running coupling, while in the large- N expansion directly in d = 3 
it is not. 

We thus consider the Lagrangian 

L = E \-^(9^ai u -d v a ¥ ) 2 + F^ (id t + a 3 a l0 + —(V-ia 3 a l ) 2 + c)F l 



1=1,2 



- g (O rl + O r2 )-g'(O p + Ol), (9) 

where I is the layer index, e 2 is the gauge coupling (which is not renormalized by 
the gauge interaction |§), Fi = (fiA, fm) T , 03 is the third Pauli matrix, and ( is the 
chemical potential. We follow the evolution of three couplings: g, g' and a = e 2 t>i?/47r 
where vf is the Fermi velocity. The one- loop renormalization of these couplings in 
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d — > 3 dimensions is 



g(fj) = 9b+ 6 2 9ba t b /3 - Cb(g 2 b + g'b) in - (10) 

= 9b + 2^/3 ~ 2c b9bg' b In - (11) 



a: 



a(/x) = a 6 (12) 

where we have used the sharp infrared cutoff at frequency \i\ the subscript b denotes 
the bare parameters; e = 3 — d, Cb = pj?/27r 2 VFb, and we do not write terms suppressed 
by powers of the ultraviolet cutoff A. The second terms on the right-hand sides 
of (0), flHD come from the gauge interaction in the BCS channel, while the last 
terms from the BCS diagrams with the local four-fermion vertices. To reiterate the 
important point, because in our model and fs have opposite gauge charges, the 
renormalization of g, g' due to the gauge interaction is of the opposite sign compared 



to the model of refs.[]T] -|Q and makes small g, g' grow in the infrared. Egs-fllTf), (11 
do not contain corrections due to the renormalization of the fields /. These are of 
the order e compared to the second terms on the right-hand sides of (PH|), ( |TTj ) and 
therefore are subleading in the (3 — d) expansion. 

We have found that the one-loop correction to four-fermion interactions from the 
gauge interaction in the 2pp channel does not contain an infrared divergence in our 
perturbation theory at d — > 3. Due to the renormalization of the fields, the particle- 
hole operators are then irrelevant as long as the bare values of couplings with 
which they are added to the Lagrangian (|^) are sufficiently small. 

Eqs.(|l0|)-(|12D also show that the naive infrared dimensions of g, g' are zero, while 
that of a is e/3 (not e, its mass dimension, as stated in ref.Qj). So, we introduce 

It is also convenient to introduce u = eg and v! = eg' instead of g and g'. Notice 



that c contains renormalizable quantity Vp but the difference between c(/i) and c& is 
inessential to the leading order of the e expansion. The one-loop beta-functions are 

du 2uX 2 /2 



+ u 2 + u' z (14) 
a in /i e7T 

cfot' 2w'A , ,„ . 

— = +2uu' 15 

a In fi err 

dX A 2 

c?ln/i 3 37r 

Note that the first two of these are singular at d — > 3 (renormalization directly at 

(i = 3 would produce In 2 fi). 

In the infrared, A runs to the stable fixed point @, A = err, which describes 

the same non- Fermi-liquid behavior as the modified fermion propagator of refs. |l], |3| 

and, probably, the fixed point of ref. 0. Substituting the fixed point value of A into 

Eqs.(|14|),(|T5|), we obtain the equations for infrared flows, t = — In/i, 

du/dt = 2u-u 2 -u' 2 (17) 
du'/dt = 2u' - 2uu' . (18) 

The resulting flows (which can be found analytically) are plotted in Fig.l. The flows 
are symmetric under u' —> —u', so we assume u' non-negative. There is a critical line 
at u = u' . If Ub > u' b , the couplings run to the infrared stable fixed point u = 2, 
u' — 0. This represents a non- Fermi-liquid normal state with repulsive in-layer four- 
fermion interaction. The interlayer electric charge transport in this state, either by 
single or pair tunneling, is suppressed. If Ub < u' b , the couplings run without limit, 
so we expect that the system generates a mass gap. Because the couplings run to 
negative values of u (attraction), we expect that this phase is superconducting. It 
should be similar in properties to that of ref.|l3|]. Because the mass gap cuts off the 
infrared renormalization, the single-particle and the particle-hole tunneling operators 
are partially recovered in this state. 



The microscopic mechanism of superconductivity may be purely electronic (in 
which case Ub is the Coulomb pseudopotential usually known as /i), or a value of Ub 
smaller than u' b may be achieved by a partial compensation of the Coulomb repulsion 
by some attractive effect, due to phonons, spin fluctuations etc. 

By analogy with the BCS theory, we expect that the critical line (in the BCS 
case, point) moves with temperature so that more points in the (u, u') plane fall 
in the domain of attraction of the normal phase. For a given point (ub,u' b ), the 
temperature at which the line passes through that point is the critical temperature 
of the superconducting to normal transition. 

Because of the singularity at e — > in the beta-functions (|14|) , (|i~5|) , the first, 
linear in u, u' terms in ([17]), (HI) are not proportional to e, and the fixed point value 
of u is not suppressed by e. This does not mean a breakdown of the e expansion. All 
the terms on the right-hand sides of Eqs.(|14D, (|15|) arise from the exceptional, BCS 
ladder diagrams, which only contribute to the beta-functions at one loop (cf. a similar 
discussion of the 1/N expansion in ref.|[].) In particular, the usual argument from 
the Fermi-liquid theory shows that there will be no terms of higher order in u, u' not 
accompanied by powers of e 2 in any of the three beta-functions. The renormalization 
of A is restricted even further. The renormalization of vp is determined from the 
fermion self-energy. Because in the (3 — d) expansion the gauge interaction adds only 
logarithmic singularities, it dressing vertices and internal lines of infrared-convergent 
Fermi-liquid graphs cannot turn them into divergent ones. Corrections to the gauge 
vertex are suppressed at low energies in d > 2 ||. Then, the only self-energy graphs 
that produce an infinite renormalization of Vp are those in which a gauge line encloses 
the rest of the graph. These are summed up by the one- loop beta-function (|16|); the 
finite renormalizations replace vpb with some vp in the bare coupling By appealing 
again to the Fermi-liquid theory, we argue that fermions with local interactions can 
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be integrated out with essentially the same effect as those without it ||, so that 
the gauge coupling e 2 is still not renormalized. Hence, we argue that the beta- 
function ( [TED for A is in fact exact. The fact that the fixed point value of u does 
not vanish at d — > 3 suggests the existence of a non-Fermi-liquid normal state in 
more isotropic three-dimensional metals close to an AF phase boundary, in which the 
quasiparticle weight vanishes logarithmically at small energies but the four-fermion 
repulsion remains constant. 

In summary, we have studied the effect of interlayer tunneling in the 2d gauge 
theory with two oppositely charged species of fermions, arising as a low-energy energy 
description of quasi-2d paramagnetic metals close to an AF phase boundary. In the 
leading order of the (3 — d) expansion, we have found two phases, superconducting and 
non-Fermi-liquid normal, suggesting that our theory may apply to the copper-oxide 
superconductors. 
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by the U.S. Department of Energy under grant DE-FG02-91ER40681 (Task B), by 
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FIGURE CAPTION 

FIG. 1. Renormalization-group flows at the fixed point value of A for the model in 
text. The arrows point towards the infrared. 
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Figure 1 



